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Abstract. The recently developed self consistent ab initio lattice dynamical method 
(SCAILD) has been applied to the high temperature bcc phase of La and Th which 
are dynamically unstable at low temperatures. The bcc phase of these metals is found 
to be stabilized by phonon-phonon interactions. The calculated high temperature 
phonon frequencies for La are found to be in good agreement with the corresponding 
experimental data. 
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1. Introduction 

The unique physical and chemical properties of the actinides and lanthanides have 
attracted interest for decades. Several reviews have been written on the subject (see, 
e.g. Refs.[Tl El E]), identifying different properties which make these two series of the 
Periodic Table unique. The uniqueness stems from the progressive filling of the f-shell as 
one traverses each series. In general one observes that the 4/-electrons in the lanthanides 
behave in a localized way, much in the same way they do in the atomic configurations. 
For the actinides (where the 5/ shell becomes filled) a more complex behavior is observed: 
for the latter part of the series (from Am and onward) the 5/ states are localized, whereas 
for the earlier part of the series (from Th to Np) the 5/ states are itinerant and form 
energy bands. This becomes evident when, e.g., considering the cohesive energy and 
equilibrium volumes [1] as well as the structural properties |5], [6]. The element Pu sits 
right on the border between these parts of the actinide series, demonstrating transition 
from itinerant /-electrons in a-phase to the localized one in 5-phase[71 E], and hence the 
electronic structure and many other properties of this material are extremely complex. 

An experimental fact which is common for the lanthanides and actinides is that 
these elements do not usually crystallize in the body centered cubic (bcc) structure 
at low enough temperatures, but they melt frequently out of this structure. In this 
paper we have chosen to investigate this complex behavior, by studies of representative 
materials of the two series, i.e. La and Th. At low temperatures La and Th have 
the double hexagonal closed packed (dhcp) and face centered cubic (fee) structure, 
respectively. The low-temperature structural properties are perfectly understood from 
first principles theory [HI [10], showing that it is the filling of the d-band that is responsible 
for the structure of La[9], whereas Th is an intermediate system which has been shown 
to have its crystal structure governed by a balance of the occupation of the 5/ and the 6d 
statesfTT]. Hence the two elements chosen in this study represent different structural and 
electronic properties in the ground state phases. Despite these marked differences these 
materials exhibit a typical behavior of the 4/ and 5/ series, namely that they melt out 
of the bcc crystal structure. It is the purpose of this report to show that the interaction 
between different phonon modes is responsible for the dynamical stabilization of the 
high temperature phase of these two series. 

2. Details of the calculations 

The lattice dynamical calculations of La and Th at elevated temperatures were 
performed using the recently developed self-consistent ab initio lattice dynamical method 
(SCAILD)|12]. In the SCAILD method one displaces the atoms of a supercell along 
several phonon modes simultaneously, in order to include the effects of phonon-phonon 
interactions upon the phonon frequencies at elevated temperatures. The simultaneous 
presence of several frozen phonons in the supercell introduces geometric disorder, i.e., 
entropy, which in turn re-normalizes the phonon frequencies. In the SCAILD approach. 



Dynamical stabilization 



3 



the atomic displacements Ur of the atoms located at equilibrium Bravais lattice sites 
R, are taken as superpositions of all the phonon modes s, with wave vectors q found to 
be commensurate with the cell, i.e, 

Here are the eigenvectors of the dynamical matrix 

R 

where $fR) are the force constant matrices and 
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where n{x) = l/(e^ — 1) is the Planck function, M is the mass of atoms, and Qq^ are the 
canonical phonon operators appearing together with the canonical phonon momentum, 
Pqs, in the harmonic Hamiltonian, Hh = Xlq,s K^q«^-q« + ^qs2q^2-qJ- 

The SCAILD method alternates between setting up atomic displacements based 
on phonon frequencies and evaluating the phonon frequencies from ab initio calculated 
forces acting on the displaced atoms. For the first iterative step, the forces stem from 
a direct force method calculation, see, e.g., Ref.[T3l HI]. The phonon frequencies and 
eigen vectors corresponding to commensurate wave vectors q, serve to calculate a set 
of atomic displacements Ur through Eqn.Q and A first principles calculation 
provides the Hellman-Feynman forces acting on the displaced atoms, and a new set of 
phonon frequencies are obtained from the Fourier transform Fq of the forces, 
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To clarify how the interactions between different phonon modes affect the 
renormalization of the phonons, we here also give the alternative expression for the 
new set of phonon frequencies ci)qs in terms of anharmonic potential coefficients 
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where <I>q,/3^(R, Ri, R2) are the third order coefficients appearing in the anharmonic part 
of the lattice dynamical hamiltonian, see for instance Ref. [IHl [16]. These coefficients 
describe, in terms of anharmonic changes to the nuclear potential energy, the underlying 
electronic structure response to the atomic displacements. Even though only third 
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order anharmonic terms are explicitly written out in Eqn. (|5]), it is now apparent 
that the simultaneous excitations of several commensurate phonons in the SCAILD 
calculation, renormalizes the phonon frequencies through anharmonic terms such as 
"^(Q) Qi) Q2, s, si, S2). Thus these anharmonic terms can be viewed as the mediators of 
the phonon-phonon interaction. 

The symmetry of the frequencies obtained through Eqn. ^ is restored through 

= ^ E (7) 

where S{q) is the symmetry group of the wave vector q and rrig the number of elements 
of the group. From these symmetry restored frequencies a new set of frequencies 

1 

<iN,) = —J2^l{^), (8) 

^ i=i 

are calculated. Here Qqsi'i), i = are the symmetry restored frequencies from 

the Nj previous iterations. The frequencies obtained through Eqn. ([s]) are then used 
to calculate a new set of atomic displacements Ur, through Eqn. ([T]) and (|3|, which in 
turn serve to calculate a new set of forces. The iterative loop is terminated when the 
frequencies obtained through Eqn. ([s]) have converged. For a more detailed discussion 
of the SCAILD method we refer the reader to the work presented in Ref. [161 E]- 

From the converged commensurate phonon frequencies and their corresponding 
eigen vectors, the force constant matrices $(R) of La and Th were calculated by inverse 
Fourier transformation. The force constant matrices were then used to calculate the 
dynamical matrices T>{q) on a mesh of 22776 q-points in the irreducible Brillouin zone 
(100x100x100 Monkhorst-Pack mesh). By diagonalization of the dynamical matrices on 
this mesh and by applying a gaussian smearing of 0.05 THz to each phonon frequency, 
the phonon density of states (DOS) was obtained for bcc-La and bcc-Th. 

Regarding the details of the force calculation we used the VASP package [TH] , within 
the generalized gradient approximation (GGA). The PAW potentials used required 
energy cutoffs of 151 eV for La and 200 eV for Th. Monkhorst-Pack grids with a 6x6x6 
mesh were used together with 0.2 eV Methfessel-Paxton smearing in all the calculations. 
The supercells used were all 64 atom cells, obtained by increasing the bcc primitive cell 
4 times along the 3 primitive lattice vectors. 

In order to get a quantitative measure of the quality of the PAW potentials used, 
the equilibrium lattice constants and bulk moduli were calculated and compared to 
their respective experimental values. Table [T] shows the results of these calculations 
together with experimental data. From these results no immediate doubt can be raised 
concerning the electronic structure part of the La and Th calculations, since calculated 
equilibrium properties are in a good agreement with experimental observations. 
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Figure 1. The right hand panel, finite temperature SCAILD calculation of La and Th 
at 1163 K and 1636 K respectively. The left hand panel, standard direct force method 
calculation [T3] of La and Th. Here the dashed curves in the left hand panel indicate 
imaginary frequencies. The filled circles and error bars are the experimental data for 
La of Ref. [TH] measured at 1163 K. 
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Table 1. The zero temperature theoretical equilibrium volumes Vq (A^), bulk moduli 
B (GPa) and pressure derivatives of the bulk moduli B' = dB/dP of the metals La and 
Th, here presented together with experimental room temperature data. The structures 
corresponding to the experimental data are hep and fee for La and Th, respectevely. 
"^Ref. [21 , ''Ref. [22 , =Ref. [23] 
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Figure 2. (color onlinc)Thc phonon density of states for bcc La and bee Th. The black 
curves are the 1163 K and 1636 K SCAILD calculations for La and Th, respectively. 
The dashed red line is the experimental bcc La data calculated from the 1163 K 
measurements of GiithofF et al [12] • Both the experimental and the theoretical density 
of states g{u!) have been normalized so that g{uj)<kj — 1. 



3. Results 

Figure [T] shows the calculated phonon dispersions for the bcc phases of La and Th at 
temperatures 1163 K and 1636 K, respectively. In the case of La, the corresponding 
experimental data of Stassis et al [.19] is displayed. Furthermore, we also observe that 
the present zero temperature calculations of La, reveal the same dynamical instabilities 
as in the previous work of Persson et al [20] . The finite temperature calculations predict 
the stability of the bcc phase of both La and Th by promoting the frequencies of the 
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phonons along the F to and T to H symmetry lines and around the P symmetry point, 
from imaginary to real. The finite temperature calculations of the bcc-La phonons result 
in an overall quantitative agreement with experimental values. Smaller deviations are 
observed around the H, P and N point of the Brillouin-zone, most likely due to finite 
size effects of the supercell used in the calculations. The frequencies at commensurate 
wave vectors are affected by the limited cell size due to the dependence on the number 
of interacting phonons (the calculation only takes into account interactions between 
phonons with commensurate wave vectors). This differs from the direct force method, 
where the cell size affects the interpolation of phonon frequencies between commensurate 
wave vectors, but not the frequencies at commensurate wave vectors. For an example of 
size effects in SCAILD calculations see Ref. dl]. From Eqn. Q and Q it can be seen 
that, at high enough temperatures, the relative shift of the squared frequency is to first 
order inversely proportional to the harmonic frequency: Su'^^/u'^g ~ a;~^^. Hence the 
transverse phonon modes, with their generally smaller phonon frequencies (compared to 
the longitudinal modes) are more strongly affected by the phonon-phononon interaction. 

Figure |2] shows the calculated DOS of bcc-La and bcc-Th at 1163 K and 1636 K, 
respectively, together with the experimental 1163 K data for La of Giithoff et al P^ . 
Here the most obvious discrepancy is that the experimental spectrum is ~ 0.5 THz 
broader than the theoretical one. Also, compared with the theoretical DOS of La the 
corresponding experimental DOS has considerable more weight located at ~ 0.8 THz. 
This lack of agreement between theory and experiment is however not surprising, since 
the electronic part of the bcc-La calculation has been performed without spin-orbit 
coupling (only scalar relativistic). 

Application of the SCAILD method in the calculation to the lattice dynamical 
properties of the rare earths La and Th, shows that the bcc structure of these elements 
is dynamically stabilized by phonon-phonon interactions. 
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